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In a recent Letter [1], Mishchenko has presented a cal- 
, culation of the minimal conductivity of the clean graphene 
' taking into account the l/r Coulomb interactions to first 
, order in perturbation theory, and questioned the results 
' obtained earlier by us [2] . He performed this calculation 
\ in three different ways and obtained three different results 
' unless, as he claimed, a proper short distance cutoff proce- 
, dure restricting large momentum transfers is implemented. 

In this comment we wish to clarify and correct some of 
' the statements made in [1] . 

First, contrary to the claims in [1] the conductivity cal- 
\ culated using the Kinetic equation is identical to the con- 
ductivity calculated from the Kubo formula. When calcu- 
, lating transport properties, it is crucial to ensure that the 
' conservation laws are built into the approximation used for 
\ the vertex. For the self-energy calculated at the Hartree- 
' Fock level, the four particle correlators must be calculated 
\ by summing the RPA bubble and ladder diagrams (see 
Fig. 1 in [1]). Such non-perturbative treatment, while 
, strictly speaking uncontrolled, does reproduce the leading 
' order perturbative expansion and does not violate Ward 
identities associated with the standard conservation laws. 
In general, however, it does violate RG scaling laws. To 
calculate the electrical conductivity, it is useful to define 
(in the reciprocal space) 

where r2„ — 2mT/ (3 is a bosonic Matsubara frequency and 
(3 = 1/T, and the vertex matrix in the real space (and 
imaginary time) is 



Setting e = vp = h = \, the Dirac particle electrical 
current operator is J^{y) = 4)'^ {y)(Jfj_\p{y). Within the pre- 
viously mentioned approximation we have 
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The above corresponds to 4 coupled integral equations. A 
rather straightforward calculation shows that the vertex 
matrix function A^(k, iVLn) can be written in the following 
form 

A^(k; iVl) = {5^,^ - k^kj^ cr^f]^{in) + ef_,^k^(73g]^{in). (4) 

Here, fi runs from 1 to 2 and the sum over repeated indices 
is understood; 6^,^ is the completely antisymmetric rank 
2 tensor and fc is a unit length wavevector. This decom- 
position allows us to reduce the set of 4 coupled integral 
equations to a set of 2 coupled integral equations for / 
and g with a non-separable vertex. At T = 0, this set of 
equations reads 
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where the Hartree-Fock energy is given by 
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and (restoring the physical units) the Coulomb potential 
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Eqs. (5) constitutes the leading approximation to the 
quantum transport theory for massless Dirac fermions in 
the collisionless limit. The electrical conductivity^ for TV 
two components Dirac flavors, in units of Ne^/h, can be 
calculated using 



cr = —2tt lim — 



d^k 
(2^)2 
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The Figure shows the numerical result found perturba- 
tively for small a — /{ehvp) which we compare with the 
analytic result found to first order in a. 

Second, and more importantly, while the conductivity 
calculated using Kubo formula would superficially indeed 
seem to depend on the cutoff procedure, further analysis 
shows, of course, that it does not. This point can be amply 
illustrated by using two different cutoffs: one hard cutoff 
A which restricts the fermion modes to the vicinity of the 
Dirac point, and second M which restricts the momentum 
transfer, i.e., 
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In the high energy parlance the second way of cutting off 
the theory corresponds to the Pauli-Villars regularization. 
The result obtained using this procedure is the dc conduc- 
tivity 



TT TT / 25 

— h a — TT 

16 6 



1 



1 



M2 



(10) 



where a = jlivp and = 4 is the number of fermion 
flavors in graphene. Naively, this would seem to lead to a 
one parameter family of results that interpolate between 
our original result (recovered when M/A — > co) and the 
one advocated by Mishchenko (when M/A 0). The root 
of this apparent ambiguity may be easily understood by 
considering the limit M ^ A, for example. One finds. 
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Fig. 1: (Upper panel) This is the result of the numerical solu- 
tion of the integral equation (5), integrated over dkk for dif- 
ferent values of e^, (vf ~ l,e ~ 1), and (Euclidean) fi. The 
analytical result for the = 0, (| — ^ tan~^ n) ^1^° plotted 
with a solid line. The agreement at = is to the 6t/i sig- 
nificant digit. The f2-slope of these lines is found at difi'erent 
values of f2 and plotted vs. (not shown). (Lower panel) The 
a-slope is then plotted for different values of fi. The red dot is 
the value found analytically. The deviation at small frequency 
is likely a result of numerical errors. 



to the Gaussian conductivity in the same system of Dirac 
fermions interacting via short-range interactions only, the 
result would be 



The last term proportional to aK/M thus corresponds to 
the perturbatively irrelevant short range part of the in- 
teraction. Indeed, if one would compute the correction 



^The conductivity can be calculated without analytical continua- 
tion by subtracting the zero frequency component of / before taking 
the limit SI — » 0: 
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where a, 6, and c are numerical constants. A is the dimen- 
sionfuU short range interaction coupling constant with the 
dimension one, that can be identified with a/M in the 
preceding discussion. The first two terms explicitly con- 
tain the cutoff, and thus clearly represent artifacts of the 
calculation that violates gauge invariance. Both of these 



terms must be dropped, so that the final result becomes 

= [cXQ + 0{fl^,X^)]j. 

This then agrees with the irrelevance of the short range 
interaction by power counting. The same result is ob- 
tained using the dimensional regularization of Veltman 
and t'Hooft, which was designed precisely to automati- 
cally discard the spurious terms of the type of the above. 

The very last term in Eq. (9) therefore is the conse- 
quence of the incorrect inclusion of the irrelevant short 
range interaction into the leading correction to the Gaus- 
sian result. Such a leading correction can come only from 
the single marginally irrelevant coupling in the theory, 
which is the Coulomb interaction. Short range interaction 
provides only the further, sub-leading corrections, propor- 
tional to the frequency. Kubo formula this way leads to 
the unique result obtained previously by us [2]. 

The last point raised by Ref. [1] is that the conductivity 
calculated using density-density correlation function fol- 
lowed by the continuity equation is divergence free and 
coincides with the Pauli-Villars result. While there is in- 
deed no logarithmic divergence in vertex correction dia- 
gram, the self energy is logarithmically divergent, as it 
may be readily seen from Eq. (12) in Ref. [1], and must be 
cut off. We admit that the agreement between these two 
procedures is puzzling. Furthermore, the Coulomb correc- 
tion to the dc conductivity calculated using the density- 
density correlator is not unique, i.e., it depends on the 
procedure implemented to regularize integrals. Namely, 
the self-energy diagram when calculated within the di- 
mensional regularization scheme, together with the ver- 
tex diagram, yields the correction to the conductivity 
CTp™"'^**^ = (7oa(ll — 37r)/6, different than one obtained 
in [1] using sharp momentum cutoff. Here ao = 7r/2 is 
the noninteracting (Gaussian) conductivity in units e^//i. 
The origin of this non-uniqueness is unclear at the mo- 
ment, but we suspect that it may be the non-gauge invari- 
ant contribution to the conductivity which is not properly 
treated within the density polarization approach. On the 
other hand, as we have already discussed, the Kubo for- 
mula yields a unique, regularization independent Coulomb 
correction to the dc conductivity found by us in Ref. [2] . 

REFERENCES 

[1] Mishchenko E. G., Europhys. Lett. 83, (2008) 17005. 
[2] Herbut I. F., Juricic V. and Vafek O., Phys. Rev. Lett. 
100 (2008) 046403. 



